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Abstract 

A  new  type  of  field-cell  interaction,  “Intracellular  Electromanipulation”,  by  means  of 
nonthermal,  wideband  electromagnetic  radiation  is  being  studied.  It  is  based  on 
capacitive  coupling  to  cell  substructures,  and  has  therefore  the  potential  to  affect 
transport  processes  across  subcellular  membranes.  This  was  verified  experimentally  by 
applying  electrical  pulses  of  60  nanosecond  duration,  with  electric  field  amplitudes  of 
>50  kV/cm  to  human  eosinophils  in  vitro.  Besides  causing  poration  of  intracellular 
membranes  without  disrupting  the  outer  cell  membrane,  sub-microsecond  pulses  were 
found  to  induce  apoptosis  in  cells.  This  effect  is  being  explored  for  cancer  treatment. 
Experiments  with  sub-microsecond  electric  field  exposure  on  a  mouse  model  showed 
reduced  tumor  growth.  In  addition  to  wideband  radiation  effects,  the  response  of  cells  to 
narrowband  electromagnetic  radiation  in  selected  millimeter-wave  (41  GHz  -  43  GHz), 
and  radio-wave  (800  kHz  to  850  KHz)  ranges  was  studied  at  moderate  power  levels 
(mW/cm2  to  W/cm2).  The  results  of  experiments  in  both  wavelength  ranges  did  not 
indicate  any  resonance  effects.  A  third  field  of  bioelectric  research,  bacterial 
decontamination  using  nanosecond  pulsed  elctric  fields  and  pulsed  corona  discharges  in 
water,  has  been  initiated.  First  results  indicate  increased  decontamination  efficiency  with 
cold  plasmas  in  water  compared  to  pulsed  electric  fields. 
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Summary 

The  program  in  Basic  and  Applied  Bioelectrics  at  Old  Dominion  University  and  Eastern 
Virginia  Medical  School,  which  was  initiated  and  supported  by  this  grant  contains  a 
medical/  biological  and  an  environmental  component.  Research  efforts  have 

a)  focused  on  the  effect  of  high  power,  ultra-fast  electrical  pulses  (wideband 
radiation)  and  that  of  moderate  power,  cw  narrowband  radiation  on  mammalian 
cells  and  tissues. 

b)  Experiments  on  the  effect  of  pulsed  electric  fields  on  bacteria  either  by  direct 
irradiation  or  through  electrical  discharges  have  been  performed. 

The  bioelectric  techniques  which  are  based  on  the  ultrafast  (nanosecond)  pulse-cell 
interactions  are  based  on  a  predicted,  and  recently  observed  effect  that  ultrashort  pulses 
(wideband  electromagnetic  radiation)  permeate  the  outer  membrane  and  affect  cell 
substructures.  One  of  the  consequences  of  this  effect  is  the  electrically  triggered 
programmed  cell  death:  apoptosis.  The  research  in  this  field  has  stimulated  a  strong 
interest  in  electro-bio-technology,  as  indicated  by  international  attendance  at  the  second 
conference  on  this  topic,  ElectroMed2001,  which  was  chaired  by  one  of  the  Co-PIs  on 
this  grant,  Stephen  Beebe.  Other  indicators  for  the  interest  of  the  scientific  community  in 
Bioelectrics  is  the  large  number  of  invited  talks,  among  them  a  plenary  talk  at  a  major 
conference  (PPPS2001).  The  enclosed  review  paper  based  on  this  talk  contains  a 
summary  of  the  status  in  modeling  and  experimental  research  in  this  new  field  of  pulsed 
electric  field-cell  interactions. 

Besides  using  wideband  radiation  to  affect  cell  structures  and  functions,  the  response  of 
narrowband  electromagnetic  radiation  on  cells  was  explored.  Since  resonance-like 
biological  effects  of  millimeter  wave  radiation  at  frequencies  of  approximately  42  GHz 
on  the  growth  rate  of  E.  coli  and  on  DNA  had  been  reported  in  the  scientific  literature,  we 
have  focused  on  this  frequency  range.  The  growth  rate  and  the  absorptior  spectrum  of  E. 
coli,  irradiated  by  millimeter  waves  in  the  frequency  range  from  41  GHz  to  43  GHz  was 
measured,  and  the  effect  of  this  radiation  on  DNA  was  studied  by  measuring  plasmid 
transformation  efficiency.  Both  the  growth  rate  variations  with  varying  frequency  and  the 
variations  in  the  result  of  the  plasmid  transformation  efficiency  experiments  were  found 
to  be  statistically  insignificant.  Resonance-like  absorption  features  observed  in  the 
absorption  spectrum  of  E.  coli  were  identified  as  modes  generated  in  the  millimeter  wave 
system,  when  the  sample  was  inserted.  The  experimental  results  indicate  that  resonance 
effects  are  unlikely  in  this  particular  frequency  range.  The  millimeter  wave  study  was 
complemented  by  an  experimental  study  of  radio-frequency  effects  on  gene  expression. 
Effects  on  gene  expression  were  determined  by  transient  transfections  of  human  cells 
with  two  different  reporter  plasmids.  The  results  of  this  study  showed  that  there  were  no 
effects  in  gene  expression  after  exposure  to  800,  835,  and  847.74  MHz  at  electric  field 
strengths  of  0.83, 1.67,  and  2.64  V/cm. 

A  third  field  of  research  supported  by  this  grant  is  bacterial  decontamination  of  liquids. 
We  have  focused  on  the  effect  of  ultrawideband  pulses  on  bacteria.  These  pulsed 
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electric  field  studies  have  recently  been  extended  to  cold  plasma  decontamination 
methods,  where  the  plasma  is  generated  in  the  liquid  itself.  Results  with  pulsed  corona 
discharges  in  water  have  indicated  a  high  efficiency  of  bacterial  decontamination  by 
means  of  cold  plasma  in  liquids.  The  leading  role  of  the  research  team  at  Old  Dominion 
University  in  pulsed  field  bacterial  decontamination  is  demonstrated  by  the  invitation  to 
write  a  review  paper  on  this  topic  in  the  IEEE  Transaction  on  Dielectrics  and  Insulation, 
and  the  appointment  of  the  P.I.  of  this  project  as  Chair  of  the  Biodielectrics  Committee  of 
the  IEEE  Dielectrics  and  Electrical  Insulation  Society. 
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Achievements  (manuscripts  related  to  the  particular  achievements  are  listed  in  italics) 


•  overview  on  status  of  “bioelectrics”  research 

Karl  H.  Schoenbach,  Sunao  Katsuki,  Robert  H.  Stark,  Stephen  Beebe,  and 
Stephen  Buescher ,  “ Bioelectrics  -  New  Applications  for  Pulsed  Power  Technology,  ” 
to  appear  in  IEEE  Trans.  Plasma  Science. 

Karl  H.  Schoenbach,  Robert  H.  Stark,  Stephen  Beebe,  and  Stephen  Buescher, 
Bioelectrics  -  New  Applications  for  Pulsed  Power  Technology,  ”  Proc.  PPPS2001, 
Las  Vegas,  NV,  June  2001. 


•  concept  of  intracellular  poration,  and  first  demonstration  of  intracellular 
poration  without  permanent  damage  to  outer  cell  membrane 

Karl  H.  Schoenbach,  “ Deactivation  of  Microbial  Pathogens  with  Intense, 
Ultrashort  Electrical  Pulses  and  High  Power  Microwave  Radiation,  ”  Proc.  1998 
ERDEC  Scientific  Conf.  on  Chemical  and  Biological  Defense  Research,  Aberdeen 
Proving  Ground,  MD,  November  1998,  p.  583. 

Karl  H.  Schoenbach,  Stephen  J.  Beebe,  and  E.  Stephen  Buescher,  “ Intracellular 
Effect  of  Ultrashort  Pulses,  ”  J.  Bioelectromagnetics  22,  440  (2001). 


•  induction  of  apoptosis  (programmed  cell  death)  in  mammalian  cells 

S.J.  Beebe,  P.M.  Fox,  LJ.  Rec,  K.  Somers,  R.H.  Stark  and  K.H.  Schoenbach 
“ Nanosecond  Pulsed  Electric  Field  (nspef)  Effects  on  Cells  and  Tissues:  Apoptosis 
Induction  and  Tumor  Growth  Inhibition,  ”  to  appear  in  IEEE  Trans.  Plasma  Science. 

S.J.  Beebe,  P.M.  Fox,  LJ.  Rec,  K.  Somers,  R.H.  Stark  and  K.H.  Schoenbach 
“ Nanosecond  Pulsed  Electric  Field  (nspef)  Effects  on  Cells  and  Tissues:  Apoptosis 
Induction  and  Tumor  Growth  Inhibition,  ”  Proc.  PPPS2001,  Las  Vegas,  NV,  June 
2001 

J.  Deng,  K.H.  Schoenbach,  E.S.  Buescher,  and  S.J.  Beebe,  “ From  Electroporation 
to  Apoptosis:  the  Effect  of  Intense,  Submicrosecond  Electrical  Pulses  on  Cells,  ” 
submitted  to  J.  Biophysics 


•  studies  on  the  effect  of  microwave  (41  GHz  to  43  GHz)  and  radiowave  (800  kHz 
to  850  kHz)  irradiation  on  E.  coli,  DNA,  and  mammalian  cells 

G.  Yu,  E.  Coin,  K.H.  Schoenbach,  M.  Gellerman,  P.  Fox,  L.  Rec,  S.  Beebe,  and  S. 
Liu,  “A  Study  on  Biological  Effects  of  Low-Intensity  Millimeter  Waves,  ”  submitted  to 
IEEE  Trans.  Plasma  Science. 

Elizabeth  A.  Coin  “A  Study  of  Radiofrequency  Effects  on  Gene  Expression”,  M.S. 
Thesis,  December  10,  2001 (Old  Dominion  University,  Department  of  Electrical  and 
Computer  Engineering) 
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•  modeling  of  the  pore  kinetics  in  cell  membranes  under  the  influence  of 
ultrashort  electrical  pulses 

R.P.  Joshi  and  K.H.  Schoenbach,  “ Electroporation  Dynamics  in  Biological  Cells 
Subjected  to  Ultrafast  Electrical  Pulses:  A  Numerical  Simulation  Study,  ”  Phys.  Rev . 
E  62, 1025  (2000). 

R.P.  Joshi,  Q.  Hu,  R.  Aly,  K.H.  Schoenbach,  and  H.P.  Hjalmarson,  “Self- 
Consistent  Simulations  of  Electroporation  Dynamics  in  Biological  Cells  Subjected  to 
Ultrashort  Electrical  Pulses, ”  Phys.  Rev.  E  64,  (2001) 

R.P.  Joshi,  Q.  Hu,  K.H.  Schoenbach, .  and  H.P.  Hjalmarson,  “Theoretical 
Predictions  of  Electro-Mechanical  Deformation  of  Cells  Subjected  to  High  Voltages 
for  Membrane  Electroporation, "  to  appear  in  Phys.  Rev.  E. 


•  studies  on  pulsed  electric  field  effects  on  bacteria 

Karl  H.  Schoenbach,  Ravindra  P.  Joshi,  Robert  H.  Stark,  Frederick  Dobbs,  and 
Stephen  J.  Beebe,  “Bacterial  Decontamination  of  Liquids  with  Pulsed  Electric 
Fields,  ”  ReviewPgper,  IEEE  Trans,  on  Dielectrics  and  Electrical  Insulation  7,  637 
(2000). 

R.  E.  Aly,  R.  P.  Joshi,  R.  H.  Stark,  K.  H.  Schoenbac,  and  S.  Beebe,  “ The  Effect  of 
Multiple,  Microsecond  Electrical  Pulses  on  Bacteria,  ”  Proc.  PPPS2001,  Las  Vegas, 
NV,  June  2001. 


•  bacterial  decontamination  of  water  using  cold  plasmas 

A.  Abou-Ghazala,  S.  Katsuki,  K.H.  Schoenbach,  F.C.  Dobbs,  and  K.R.  Moreira, 
“Bacterial  Decontamination  of  Water  by  Means  of  Pulsed  Corona  Discharges, " 
submitted  to  IEEE  Trans.  Plasma  Science. 

A.  Abou-Ghazala,  S.  Katsuki,  K.  H.  Schoenbach,  F.  C.  Dobbs,  and  K.  R.  Moreira, 
“ Bacterial  Decontamination  of  Water  by  Means  of  Pulsed  Corona  Discharges,  ” 
Proc.  PPPS2001,  las  Vegas,  NV,  June  2001. 


•  development  of  pulsed  power  systems  for  bioelectrics 

K.H.  Schoenbach,  R.H.  Stark,  J.  Deng,  R.  El-Sayed  Aly,  S.J.  Beebe,  and  E.S. 
Buescher,  “Biological/Medical  Pulsed  Electric  Field  Treatments,"  Conf.  Record, 
2000  Twenty-Fourth  Intern.  Power  Modulator  Symposium,  June  2000,  Norfolk,  VA, 
p.  42. 

J.  Deng,  R.H.  Stark,  and  K.H.  Schoenbach,  “A  Nanosecond  Pulse  Generator  for 
Intracellular  Electromanipulation,”  Conf.  Record,  2000  Twenty-Fourth  Intern. 
Power  Modulator  Symposium,  June  2000,  Norfolk,  VA,  p.  47. 


•  organization  of  the  Second  International  Symposium  on  Nonthermal 
Medical/Biological  Treatments  Using  Electromagnetic  Fields  and  Ionized  Gases, 
“Elect  roMed2001”  in  Portsmouth,  VA,  May  21  -  23, 2001 

Symposium  Record  Abstracts 
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S.  Beebe:  Chairman 


•  sixteen  invited  presentation  on  the  topic  of  Bioelectrics,  two  of  them  plenary 
talks,  and  one  of  them  a  keynote  address. 

listed  in  the  Bibliography  section  of  this  report 


•  patent  on  “Bioelectrics” 

Method  and  Apparatus  for  Intracellular  Electromanipulation 
Patent  No:  6,326,177 
Patent  Date:  December  4,  2001 

Inventors:  Karl  H.  Schoenbach,  Stephen  J.  Beebe,  and  E.  Stephen  Buescher 
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The  temporal  dynamics  of  electroporation  of  cells  subjected  to  ultrashort  voltage  pulses  are 
studied  based  on  a  coupled  scheme  involving  the  Laplace,  Nemst-Plank,  and  Smoluchowski 
equations.  It  is  shown  that  a  finite  time  delay  exists  in  pore  formation,  and  leads  to  a  transient 
overshoot  of  the  transmembrane  potential  Vmem  beyond  1 .0  Volt.  Pore  resealing  is  shown  to  consist 
of  an  initial  fast  process,  a  1  O'4  second  delay,  followed  by  a  much  slower  closing  at  a  time  constant  of 
about  10'1  second.  This  establishes  a  time-window  for  effective  killing  by  a  second  pulse.  The 
results  are  amply  supported  by  our  experimental  data  for  E.-coli  cells,  and  the  time  constant  also 
matches  experiments.  An  electro-mechanical  analysis  for  analyzing  cell  shape  changes  is  also 
presented.  Our  calculations  show  that  at  large  fields,  the  spherical  cell  geometry  can  be  significantly 
modified,  and  even  ellipsoidal  forms  would  be  inappropriate  to  describe  the  deformation.  Values  of 
surface  forces  obtained  are  in  very  good  agreement  with  the  1-10  nN/m  range  reported  for  membrane 
rupture.  It  is  also  demonstrated  that,  at  least  for  the  smaller  electric  fields,  both  the  cellular  surface 
area  and  volume  change  roughly  in  a  quadratic  manner  with  electric  field.  Finally,  it  is  shown  that 
the  bending  moments  are  generally  quite  small  and  can  be  neglected  for  a  simpler  analysis. 
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I.  INTRODUCTION 

Electroporation  is  a  well  known  physical  process  in  biological  cells  [1-3].  It  involves  rapid 
structural  rearrangement  of  the  membrane,  in  response  to  an  externally  applied  electric  field.  The 
most  prominent  observable  effect  is  a  rapid  increase  in  the  electrical  conductivity  by  several  orders 
of  magnitude  [4].  This  is  ntiributed  to  the  formation  of  aqueous  pathways,  or  pores,  in  the  lipid 
bilayer  of  the  membrane.  The  opening  of  such  channels  (or  more  appropriately,  transient  aqueous 
pores)  enables  the  transport  of  ions  and  water-soluble  species  both  into  and  out  of  individual  cells. 
Electroporation  can,  therefore,  be  used  to  initiate  large  molecular  fluxes  for  purposes  of  introducing 
genetic  material  into  cells,  manipulation  of  cells  and  tissues,  and  other  applications  in  biotechnology 
[5-9]. 

Electroporation  has  also  been  linked  to  the  non-thermal  killing  of  micro-organisms  subjected 
to  strong  electric  fields  [10].  For  this  reason,  it  offers  great  potential  for  de-contamination  and  the 
elimination  of  harmful  micro-organisms  and  bio-hazards.  Traditionally,  most  electroporation  studies 
have  focused  on  relatively  low  external  voltages  applied  over  extended  time  periods  ranging  from 
several  tens  of  microseconds  to  milliseconds  [1 1].  Recently,  work  has  focused  on  the  use  of  much 
shorter,  high-voltage  pulses  with  electric  fields  as  high  as  100  kV/cm  [12-15]  and  durations  well 
below  the  micro-second  range.  There  appear  to  be  several  fundamental  advantages  in  using  short 
electric  pulses  for  cellular  manipulation.  First,  negligible  thermal  heating  of  the  biological  matter 
can  be  expected  to  occur  due  to  the  short  time  duration.  Much  lower  energies  are  required  for  pulsed 
inputs,  and  yet  large  values  of  the  electric  fields  and  peak  powers  can  be  obtained  [16].  Also,  pulsed 
fields  afford  a  way  by  which  the  time  scales  can  easily  be  manipulated. 

Electroporation  dynamics  were  first  treated  based  on  the  Smoluchowski  equation  by 
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from  an  external  step  voltage  are  computed  based  on  the  current  continuity  and  Laplace  equations. 
Current  flows  are  computed  corresponding  to  the  electric  fields  at  each  time  step  from  a  continuum 
Nemst-Planck  type  model  that  includes  diffusion  [27, 28, 32].  In  keeping  with  the  literature  [33-35], 
it  is  assumed  here  that  two  types  of  pores,  hydrophilic  and  hydrophobic,  exist.  Each  of  the  two  pore 
types  were  characterized  by  a  “formation  energy”  function,  E(r)  dependent  on  the  pore  radius  V . 
The  following  pore  energy  function  was  chosen  in  keeping  with  the  published  and  accepted  model 
[4,27,36,37]. 

E(r,t)  =  2  rthr  a(«>)  [  Ii(r/ro)  /  lo(r/ro)]  -rcapV2r2  ,  and  (la) 

E(r,t)  =  2  iryr  -  [  oJr  27tr(r*)r*  dr*  ]  +  (C/r)2  -  7t  ap  V2  r2  ,  (lb) 

for  hydrophobic  and  hydrophilic  pores,  respectively.  In  the  above,  Ii  and  lo  are  the  modified  BesSel 
functions  of  the  zeroth  and  first  order,  respectively,  h  is  the  membrane  thickness,  o(«>)  is  a  constant 
equal  to  5  x  1 0"2  N  m'1 ,  while  ro  represents  a  characteristic  length  scale  over  which  the  properties  of 
water  change  between  the  interface  and  the  bulk,  taken  to  equal  1  nm.  A  simple  heuristic  model  has 
recently  been  used  [38]  for  r(r)  in  which  r(r)  =  To  [1  -  r 2/r«2].  This  is  based  on  the  idea  that  tension 
is  proportional  to  the  membrane  area,  at  least  to  first  order.  The  primary  effect  of  such  variations  in 
r,  would  be  the  creation  of  an  additional  local  minima  in  the  pore  energy  function.  From  the 
standpoint  of  electroporation,  this  means  that  instead  of  expanding  indefinitely  beyond  the  unstable 
maxima,  the  pores  ultimately  stabilize  at  a  high  radial  value.  The  variable  tension  concept  had  been 
proposed  by  Winterhalter  and  Helfrich  [39],  as  well. 

Typical  values  for  the  various  parameters  are  given  in  Table  I.  The  value  of  "C"  in  (lb)  was 
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chosen  to  be  9.67x  10‘15  Jl/4  m  in  keeping  with  the  reports  by  Neu  et  al.  [34]  as  it  yields  values  close 
to  those  measured  by  Glasser  et  al.  [27].  The  coefficient  ap  is  a  property  of  the  membrane  and  its 
aqueous  environment.  In  the  simplest  continuum  approximation  [27],  and  expressed  in  terms  of  the 
membrane  thickness  "h”  and  the  permittivities  "e*”  and  "e  m"  of  water  and  the  membrane, 
respectively,  as :  ap  =  (£w-e  m)/[2h] .  The  energy  fraction  E(r)  determines  the  “drift  flux”  for  pores 
in  r-space  and  therefore,  governs  the  growth  or  contraction  of  pores  at  any  given  radius  "r".  In 
general,  the  presence  of  a  membrane  voltage  reduces  the  maxima,  and  can  even  quell  the  energy 
barrier  completely  beyond  a  critical  voltage  value.  For  transient  voltage  pulses,  stability  would 
depend  on  the  ability  of  pores  to  drift  past  the  barrier  maximas  within  the  duration  of  the  applied 
voltage  pulse.  It  was  assumed  here  that  all  pores  are  initially  created  as  hydrophobic  (i.e.  non¬ 
conducting)  at  a  rate  S(r)  per  unit  area  of  the  membrane  at  a  rate  given  by : 

S(r)  =  {(vc  h  )  /  (kB  T)}  [dE(r)/dr]  exp[-E(r)/(kBT)]  dr  dt ,  (2) 

-  where  vc  is  an  attempt  rate  density  [36],  E(r)  the  en-.rgy  for  hydrophobic  pores,  T  the  operating 
temperature,  and  kB  the  Boltzmann  constant.  If  a  non-conducting  pore  is  created  with  a  radius  r  >  r* 
(=  0.5  nm),  it  spontaneously  changes  its  configuration  and  transforms  into  a  conducting,  hydrophilic 
pore.  All  conducting  pores  then  survive  as  long  as  their  radii  remains  larger  than  r* .  Destruction  of 
a  conducting  pore  occurs  only  if  it  drifts  or  diffuses  in  r-space  to  a  value  below  r*  .  Due  to  the 
exponential  term  in  (2),  most  pores  are  created  with  very  small  radii. 

The  Smoluchowski  equation  that  governs  the  pore  dynamics  is  given  in  terms  of  the  pore 
density  distribution  function  n(r,t)  as  : 
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dn'(r,t)/dt  +  {D/[kB  T]}  [  d{E(r)  n(r,t)}/5r  ]  -  D  [a2  n(x£)ld?  ]  =  S(r  ) , 


(3) 


where  S(r)  is  the  source  term  as  given  in  (2),  while  D  is  the  pore  diffusion  constant  given  in  Table  I. 
The  process  of  diffusion  represents  a  “random  walk"  of  the  pore  radius  in  “r-space”.  Numerical 
simulations  of  the  dynamic  pore  distribution  were  carried  out  based  on  a  time-domain,  finite- 
difference  discretization  of  the  governing  Smoluchowski  equation.  An  upperbound,  rmax ,  of 2000  A 
was  set  on  the  pore  radius,  and  this  entire  r-space  was  uniformly  divided  into  5000  segments.  A 
“reflecting  boundary”  was  assigned  at  r  =  rmax ,  which  was  implemented  by  setting  the  pore  flux  to 
zero  at  r  =  rmax  .  Mathematically,  this  amounts  to  a  Neumenn  condition:  [dn(r,t)/dr  ]|  r  =  rmax  =  - 
[dE(r,t)/dr]  [n/(ke  T)]  r  =  m»ax  •  At  the  other  end,  absorbing  boundary  conditions  were  implemented 
by  setting  n(0,t)  =  0.  The  time  step  dt  =10'1 1  second  was  chosen  as  it  is  to  be  much  smaller  than  the 
fluctuation  rate  Vd  representing  the  fastest  time  constant  in  the  system  [27]. 

The  schematic  shown  in  Fig.  1  was  used  to  represent  a  cell  suspended  in  a  medium,  and  the 
relevant  equations  applied  to  this  geometry.  The  inner  spherical  region  was  taken  to  have  a  radius  “a” 
and  conductivity  a  „ .  The  applied  electric  field  Fo(t)  was  taken  to  be  along  the  z-axis.  The  cellular 
membrane  was  assigned  a  thickness  “b  -  a”  ,  while  the  outer  region  denoting  the  suspension  was 
assigned  a  conductivity  a  out  .  Due  to  axial-symmetry,  the  potentials  satisfying  the  Laplace  equation, 
could  be  expressed  as: 


Uin  (r,t)  =  A0  (t)  P0  +  Ai  (t)  r  Pi  +  A2  (t)  r2  P2  +  ...  =  2^=o,«  Aj  (t)  r*  Pj  ,  (3a) 
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umem  (r,t)  =  £j=0,oe  [ Bj (t) r* Pj  +  Cj(t)Pj/r^  ]  ,  (3b) 

and  Uout  (r,t)  =  -  F0(t)rP,  +  Zj=0,-  Dj  (t)Pi/ii+1  ,  (3c) 

where  Ui„  (r,t),Umem  (r,t),  and  Uout  (r,t)  are  the  potentials  at  the  inner-,  the  membrane- and  outer- 
regions,  Pj  is  the  j*  order  Legendre  polynomial,  and  F0  (t)  the  externally  applied  electric  field.  Also, 
Aj  (t),  Bj  (t),  Cj(t),  and  Dj  (t)  denote  the  coefficients  of  the  Legendre  series  expansions  that  can  be 
determined  by  applying  matching  boundary  conditions  at  the  interfaces  of  the  three  regions.  Here, 
the  Laplace  instead  of  Poisson’s  equation  has  been  used  on  the  assumption  that  charge  inequalities 
arising  from  ionic  transport  can  be  ignored  on  the  short  time  scales.  Invoking  continuity  in  potential 
and  current  density  then  leads  to  the  following  expression  for  the  membrane  potential  Vmem  (t)  : 

V mem  (t)  A(t)  cos(0)  [  {b-b  /a  }  {l+2aout/aj„}/  {(b/a)  2o  out / o  jn  }  -  (b-a)  ] 

+  3  cos(0)  F0  (t)  (o  out  /a  in )  /  {(b/a)3  +  2 a  ou,  /a  in }  ,  (4a) 

where  0  is  the  angle  with  respect  to  the  z-axis  (and  hence,  the  applied  electric  field  direction), 
and  A(t)  satisfies  the  following  ordinary  differential  equation  : 

CM  [dA(t)/dt  ]  [{b-b3/a2}  {l+2aou,/ain}/  {(b/a)3  +  2oout/ain }  -  (b-a)  ]  =  -  Jmem  (t)  - 
<^in  [  {  2(a0Ut/om )  {(b/a)  - 1}  A(t)  —  3  (o0Ut/ojn )  Fo  (t)  (b/a)  }/{(b/a)  +  2oout/crin }  ] .  (4b) 

Under  steady-state,  the  transmembrane  potential  Vmem  from  the  above  equation  simplifies  to:  - 1 .5 


8 


a  E  cos(0)  for  b  ~  a,  in  accordance  with  reports  in  the  literature  [40, 41]. 

The  conduction  current  density  Jmem  (t)  needs  to  be  specified  in  order  to  solve  for  the 
potentials  in  (4).  A  one-dimensional  approximation  of  the  Nemst-Planck  expression  for  ionic  flow  is 
[27,28,32] : 

I  (t)  =  it  a  R2  N(t){exp[q  Vmem(t)  /(kBT)]  -  1  }/{oJh  exp[q  Vmem(t )  { l-x/h}/(kBT)  +  w(x)]  dx} ,  (5) 

where  R(t)  is  the  pore  radius,  a  the  conductivity  of  the  aqueous  solution  that  fills  the  pore,  w(x)  the 
energy  barrier  to  ionic  flow  through  the  pores,  “h”  is  the  membrane  thickness,  and  N(t)  the  pore 
density.  A  simple  trapezoidal  form  for  w(x)  as  given  by 

w(x)  =  qA/(kBT)  {x/hi}  for  0  <  x  <  hi  ,  (6a) 

w(x)  =  qA/(kBT)  for  hi  <x  <  h-hi  ,  (6b) 

and,  w(x)  =  qA/(kBT)  {(h-x)/h,}  for  h,  <x  <  h-h,  ,  (6c) 

was  used  for  the  barrier  energy,  and  will  be  applied  here.  In  the  above,  “A”  represents  the  peak 
barrier  height  under  zero  bias.  Values  of  o  and  A  are  known  to  be  roughly  1.3  S/m  [42]  and  2.5 
Volt  [27] ,  respectively.  The  parameter  hi  is  the  length  near  the  pore  entrance  over  which  the  barrier 
profile  would  be  changing  linearly  for  an  unbiased  cell  roughly  equals  0.15  times  the  membrane 
thickness.  In  addition,  an  ionic  component,  Iion,  which  is  orders  of  magnitude  lower  in  strength,  has 
to  be  included.  This  ionic  current  density  is  [42]:  lion  =1.9  (Vmem  +  .083). 
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It  must  be  pointed  out  that  the  above  equation  (5)  is  somewhat  inaccurate  because:  (i)  The 
barrier  peak  “A”  is  assumed  to  be  independent  of  the  pore  radius.  An  approximate  form,  which  is 
correct  for  an  infinitely  long  cylindrical  pore  geometry,  for  w(r)  has  been  derived  to  be  [43] :  w(r)  = 
5  x  10'9  /  r .  In  our  calculations,  this  w(r)  function  was  explicitly  used,  (ii)  Next,  equation  (5)  treats 
the  pore  radius  to  be  a  constant.  However,  not  only  would  the  radius  of  the  pores  change  under 
transient  conditions  upon  the  application  of  an  external  voltage  waveform,  but  also  the  pores  would 
not  all  be  identical  in  size.  The  distribution  n(r,t)  as  predicted  by  the  Smoluchowski  equation,  would 
impart  a  heterogeneous  spread  in  the  “R”  parameter,  (iii)  Finally,  N(t)  which  is  an  integral  quantity 
needs  to  be  obtained  through  a  suitable  integration  of  n(r,t)  over  r-space.  Treating  N(t)  as  a 
fundamental  independent  variable  is  incorrect. 

In  order  to  correct  the  above  shortcomings  of  equation  (5),  the  following  modified  current- 
voltage  (I-V)  characteristics  were  used  in  this  study : 

I(t)=rcooJ~  [r2  n(r,t)  {exp[qVmem(t)  /(kBT)]-l  >/{0/h  exp[qVmem(t){l-x/h}/(kBT)+w(r,x)]  dx}],  (7) 

which  includes  an  integration  in  r-space  over  the  time-dependent  distribution  n(r,t).  In  the  process, 
the  role  of  inter-cellular  variations  and  the  size  distribution  of  pores  are  both  taken  into  account.  A 
distribution  of  barrier  energies  and  their  fluctuation  from  site  to  site  is  also  automatically  included  in 
this  formulation. 

B.  Cell  Stress  and  Deformation: 

Our  basic  stress  model  is  based  on  the  classical  small  deformation  theory  of  thin,  elastic 
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shells  [30].  Since  the  thickness  of  cell  membranes  is  on  the  order  of  5  pm,  compared  to  their  radii  of 
~1  pm,  the  shell  theory  is  quite  appropriate.  We  follow  the  notation  of  Flugge  [30]  for  the  stresses 
and  moments.  Thus,  N^,  denotes  the  meridional  force  per  length ,  No  the  hoop  force  per  length,  and 
the  shear.  Also,  pr ,  p*  and  pe  are  the  externally  imposed  stresses  (which  could  include  internal 
osmotic  pressure),  while  r  is  the  distance  from  the  axis  of  rotation,  rj  the  radius  of  curvature,  and  r2 
the  distance  of  intersection  of  the  radius  of  curvature  and  the  axis  of  revolution.  In  the  present 
context,  pr ,  p,|>  and  pe  will  be  non-zero  due  to  the  Maxwell  stress  tensor  that  arises  from  the  external 
field.  From  the  geometry,  r  =  r2  sin(<|>),  while  the  elemental  distance  "ds"  along  the  meridian  is  given 
by:  ds  =  ri  d<(>,  where  <|>  is  the  angle  between  a  normal  to  the  shell  and  its  axis  of  revolution.  The 
parameter  0  denotes  the  meridional  angle.  Finally,  M$  ,  M0  and  are  the  bending  moments 
(dimensions  of  force),  while  Qe  and  Q*  transverse  forces  per  length  that  arise  from  bending  theory. 
At  equilibrium,  the  balance  of  all  forces  and  moments  yields  the  following  six  equations: 

d{r  Nv}/d<|>  +  ri  d{Ne<|,}/d0  -  ri  N0  cos(<)>)  -  r  Q*  =  -  r  T\  p*  ,  (8a) 

d{r  N$e  }/d<|>  +  r,  d{Ne}/d0  +  ri  Ne^  cos(<|>)  - rj  Qe  sin(<»  =  -  r  ri  p0  ,  (8b) 

ri  Ne  sin(<|))  +  r  N^,  +  ri  dQe/d0  +  d{r  Q0}/d<()  =  r  T\  pr  ,  (8c) 

d{r  M,|,}/d<I>  +  ri  d{Me<>}/d0  -  ri  M0  cos(<|>)  =  r  rt  ,  (8d) 

d{r  M^}/d())  +  rj  d{Me}/d0  +  ri  Me^,  cos((|))  =  r  ri  Q0  ,  (8e) 

/  r  i  -  Me<|>  /  r2  —  N^e  -  Ne$  .  (8f) 

The  current  problem  of  interest  involves  an  inherent  axial  symmetry  along  the  direction  of  the 
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applied  electric  field.  For  such  axisymmetric  cases,  the  derivatives  with  respect  to  the  angle  0  drop 
out,  while  the  shearing  forces  N*e  and  N^,  the  twisting  moments  Me*  and  ,  and  the  transverse 
shear  Qe  all  vanish.  Also,  the  load  component  pe  is  zero.  Consequently,  the  following  simpler  set  of 
equations  result: 


d{r  N<t,}/d4)  -  n  Ne  cos(<}>)  -  r  Q6  =  -  r  ri  p*  ,  (9a) 

rj  Ne  sin(<|))  +  r  N<,  +  d{r  Q*}/d<f>  =  rripr,  (9b) 

d{rHt,}/d<|)  -  ri  M0cos(<f>)  =  r  r(  Q^,  .  (9c) 

The  above  three  equations  contain  five  unknowns,  and  need  to  be  supplemented  by  the  stress-strain 
relationships.  In  the  elastic  regime,  the  stresses  can  be  related  to  the  displacements  v  (transverse) 
and  w  (normal)  in  the  following  manner  [30] : 

Hi,  =  (12K/12)  [{dv/d<|>  +  w}/rj  +v  {vcos(<|))  +  wsin(<J>)}  /r]  ,  (10a) 

Ne  =  (12  K /t2  )  [{v  cos(<|>)  +  w  sin(<]>)}  / r  +  v  {dv/d<j>  +  w}/ ri  ]  ,  (10b) 

M4,  =  (K/  ri)  £d({dw/d<|>}/ ri)/d<t>  +v  {  d[cos(<|>){dw/d<l>}/r]/d<|> }]  ,  (10c) 

Me  =  (K/rj)  [cos(<|>){dw/d<J>}/r]  +  v  d({dw/d<J)}/ri)/d<j>]  ,  (lOd) 

where  K  is  the  flexural  rigidity  (i.e.  bending  stiffness),  v  the  Poisson’s  ratio,  t  the  shell  thickness 
assumed  to  be  a  constant,  and  w  and  v  the  displacements  due  to  deformation  along  the  radial  and 
angular  directions.  Equations  (10a)-(10d)  involve  the  displacements  w  and  v  that  constitute  two 
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additional  unknowns  of  the  problem.  Thus,  the  combined  set  of  equations  [9(a)-9(c)]  and  [10(a)- 
10(d)]  yield  a  system  of  seven  equations  for  the  seven  unknowns  that  can  be  solved. 

C.  Electric  Field,  Analysis : 

The  electric  fields  that  give  rise  to  the  Maxwell  stress  tensor  can  be  computed  based  on  the  Legendre 
polynomial  model  for  the  potential  distributions.  Tedious,  but  straight  forward  calculations  yield  the 
following  expressions  for  the  electric  fields  Fr(r)  and  F$(r)  just  outside  the  membrane  (i.e.  at  r  =  b+): 


F^r=b)  =  [3  F0  +  2  (C,/b3)  {1  +  (b/a)3  T }]  cos(<j>)  =  Fr  cos(<t>) ,  (lOf) 

F^(r=b)  =  {Ci  /  b3}sin(<|>)  [1  +  (b/a)3  T  ]  =  F*  sin(<|>)  ,  (lOg) 
where,  Ci  =  -  3Fo8out  /  [{T  Cmem  /a3}  - 2emem  /b3  +  {2£oUt /b3}[l  +  (b/a)3 T  ]  ,  (lOd) 
and,  T  —  {2c  mem"^£  in  }/ {£  mem  "  Gin  }  »  (10e) 

where  e  i„,  e  mem  and  e  out  are  the  permittivities  of  the  inner,  membrane,  and  outer  regions, 
respectively.  For  e  in  =  c  mem  =  e  out ,  the  above  equations  reduce  to:  U(r)  =  -  F0  r  cos(<t>),  Fr  (r)  =  F0 
cos(<t>),  and  F$(r)  =  -  Fo  sin(^>). 

The  expressions  for  the  surface  forces  N$  (<{>)  and  Ne  (<)))  become: 

N4,  (<|>)  =  [F*  /  {r2  sin2(<|>)}]  {  of1’  ri  r2  [cos(2a* )  cos«>*  )sin(<j>*  )  +  sin(2a)  sin2(<|>* )]  d<}>*  }  ,  (11a) 
No  (4>)  =  r2  F*  cos(2a)  -  [r2  /  ri]  N*  ,  where  a*  =  tan'1  [-tan(<|>*)  F^/FJ  .  (1  lb) 


m.  RESULTS  AND  DISCUSSION 
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Self-consistent  simulations  based  on  the  coupled  Laplace-Nemst-Planck-Smoluchowski 
equations  were  carried  out  next  to  evaluate  the  temporal  response  to  ultrashort,  external  electric 
pulses.  A  10  kV/cm  rectangular  external  electric  field  pulse  with  a  4  x  10-6  duration  was  assumed. 
These  parameters  were  chosen  in  keeping  with  actual  pulsed  field  experiments  conducted  on  E.  coli 
in  our  laboratory  to  facilitate  comparisons  between  theory  and  experiment.  The  cell  radius  was 
chosen  to  be  1.0  pm  and  a  membrane  thickness  of  5  nm  which  is  roughly  characteristic  of  E.  coli. 
Fig.  2  shows  the  dynamic  evolution  of  the  pore  density.  An  initial  delay  of  about  5  ns  seen  in  Fig.  2a 
is  due  to  membrane  charging  and  for  Vmem  to  build  up  to  levels  at  which  the  pore  creation  rate 
becomes  substantial.  A  peak  value  of  about  2  x  1 014  m'2  is  reached  after  about  20  ns.  Subsequently, 
the  pore  density  shows  a  slight  monotonic  decrease  over  the  remaining  duration  of  the  4  ys  external 
pulse.  This  occurs  due  to  a  substantial  increase  in  the  membrane  conductance  and  a  consequent 
decay  in  Vmem  that  controls  the  pore  generation.  Details  of  the  time  dependent  membrane  voltage  are 
shown  in  Fig.  3.  The  voltage  exceeds  1.0  volt  at  about  15  ns,  and  reaches  a  peak  value  of  roughly 
1.2  volts.  At  this  point  the  pore  conductance  increases  to  such  a  degree  that  the  voltage  across  the 
membrane  capacitance  begins  to  fall.  The  overall  result  is  a  “voltage  overshoot”  behavior.  It  agrees 
well  with  a  previous  report  on  the  time  dependent  behavior  of  the  membrane  voltage  by  Meier  et  al. 
[29].  As  the  external  electric  field  is  turned  off  beyond  4  ys,  the  transmembrane  potential  falls 
dramatically  with  a  time  constant  in  the  sub-microsecond  range.  The  fast  drop  off  is  the  result  of  a 
large  conductance,  and  hence,  a  low  internal  “RC”  time  constant.  A  final  steady  state  value  of  about 
-80  mV,  equal  to  the  rest  potential  is  finally  reached.  The  corresponding  influence  on  the  pore 
density,  as  seen  from  Fig.  2a,  is  a  sharp  decrease  by  about  fifty  percent  following  the  turn-off. 
Beyond  this,  the  density  continues  to  decrease,  but  at  a  relatively  low  rate  until  about  0.1  ms.  This 
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implies  that  many  of  the  pores  tend  to  remain  open,  well  after  the  4  ps  pulse  is  switched  off.  Hence, 
a  second  electric  pulse  applied  within  this  duration  is  likely  to  have  a  substantial  damaging  effect. 
This  is  borne  out  in  our  experimental  measurements  on  E.  coli,  as  discussed  later.  Also,  this  0. 1  ms 
time  delay  corresponds  well  with  an  experimental  report  by  Meier  et  al.  [29] .  Beyond  0. 1  ms,  the  rate 
of  pore  reduction  increases.  The  long  time  behavior  can  best  be  gauged  from  the  semi-logarithmic 
curve  of  Fig.  2b.  It  shows  a  gradual  slowing  in  the  resealing  rate.  At  the  2.4  ms  instant,  an 
“effective  decay  time  constant”  of  4.5  x  10'3  seconds  is  computed  which  leads  to  a  pore  resealing 
time  of  about  0.12  second.  In  actuality,  though,  the  duration  would  be  even  longer  since  the  decay 
curve  of  Fig.  2b  exhibits  a  continuous  slowing-down.  In  any  case,  the  resealing  values  projected 
values  are  in  the  10'1  -  102  second  range,  in  keeping  with  experimental  reports  [44, 45]. 

The  dynamical  behavior  can  easily  be  understood  in  terms  of  the  formation  energy 
characteristic.  Initially,  the  pores  that  are  not  near  the  local  minima,  and  have  values  in  r-space  that 
are  below  the  barrier.  These  pores  tend  to  drift  and  diffuse  towards  r  0,  giving  rise  to  a  fast  decay. 
However,  this  leaves  behind  an  ever-increasing  higher  fraction  of  larger  pores  that  are  near  the  local 
minima  Those  near  the  minima  move  in  r-space  primarily  through  diffusion,  and  hence,  the  system 
takes  a  long  time  to  completely  recover  to  the  original  steady  state. 

The  effect  of  including  a  second  electric  pulse  and  the  comparative  temporal  behavior  is 
shown  in  the  curves  of  Fig.  4.  The  pore  density  evolution  for  single  and  a  dual  4  micro-second,  10 
kV/cm  rectangular  electric  pulses  is  shown.  The  delay  for  the  two-pulse  situation  was  taken  to  be 
1 .2  ms.  Despite  being  outside  the  0. 1  ms  time-window,  the  application  of  a  second  pulse  is  seen  to 
have  effect  on  the  evolution.  First,  the  pore  density  is  increased  leading  to  a  short  spike.  However, 
this  increment  is  not  as  large  as  that  produced  by  the  original  pulse,  since  the  larger  pores  remain 
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open,  giving  rise  to  a  reasonably  large  membrane  conductivity.  Consequently,  the  membrane 
potential  developed  is  not  as  large,  and  the  pore  creation  rate  not  boosted  as  much.  A  second 
observation  is  that  a  two-step  decay  is  evident  following  the  second  pulse,  as  with  the  initial  pulse. 
Finally,  and  most  important,  the  decay  time  for  the  two-pulse  situation  becomes  significantly  longer 
beyond  2  ms.  This  implies  that  applying  two  or  mo.  e  pulses  can  slow  resealing  considerably.  From  a 
practical  standpoint,  this  has  applications  for  drug  delivery  over  prolonged  durations. 

For  a  more  quantitative  evaluation,  it  is  instructive  to  examine  the  pore  distribution  functions 
at  various  times.  Shown  in  Fig.  5  are  five  distribution  functions  at  various  times  for  the  single-  and 
dual-pulse  scenarios.  The  highest  curve  for  the  dual-pulse  situation  is  at  a  time  instant  of  1 .205  ms 
and  hence,  just  after  the  completion  of  the  second  pulse.  Though  the  peak  lies  at  about  1.5  nm,  a 
fairly  broad  “tail”  at  larger  radii  is  apparent.  At  a  later  time  of  t=1.23  ms  for  the  dual-pulse  situation, 
the  peak  is  seen  to  be  reduced  appreciably  and  broadened  due  to  diffusion  in  r-space.  Finally,  at  the 
longest  time  of  2.43  ms,  the  dual-pulse  distribution  is  seen  to  have  become  quite  negligible  for  small 
pore  radii,  but  has  a  well-defined  peak  at  about  52  nm.  This  location  corresponds  to  the  potential 
minima  for  the  r«  =  65  nm  curve.  It  is  thus  evident  tint  a  well-defined  population  of  relatively  large 
pores  remains  and  is  long-lived.  For  the  single-pulse  situation,  the  snapshot  distributions  at  t  =  1.2 
ms,  and  2.4  ms  are  shown.  The  distribution  is  relatively  broad  at  1.2  ms,  and  similar  to  the  dual¬ 
pulse  case,  settles  down  to  acquire  a  local  peak  at  52  nm.  However,  the  52  nm  peak  is  smaller, 
signifying  a  greater  number  of  pores  for  the  dual-pulse  case.  Beyond  these  times,  the  resealing  can 
be  expected  to  be  slow,  and  pore  may  remain  open  for  a  very  long  time. 

Some  experimental  findings  and  measured  data  for  E.  coli  cells  subjected  to  short  electric 
pulses  are  presented.  The  results,  plotted  in  Fig.  6,  show  the  measured  viability  of  E.  coli  subjected 
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to  two  4  ps,  13  kV/cm  electric  pulses  with  variable  delays.  The  delay  between  the  first  and  second 
pulses  ranged  from  5  ps  to  50  seconds.  A  sudden  jump  in  the  viability  is  seen  in  going  from  a  0.1 
ms  delay  to  a  1 .0  ms  delay  for  the  second  pulse.  This  is  in  keeping  with  the  simulation  result  of  Fig. 
3a  that  shows  a  0.1  ms  time-window  during  which  the  pores  are  mostly  open,  making  the  system 
most  vulnerable  to  destruction  by  a  second  electric  pulse.  Second,  there  is  no  significant  change  in 
the  experimental  viability  in  going  from  a  1 .0  ms  delay  to  a  1  second  delay.  This  again  is  in  keeping 
with  the  prediction  of  a  slow  pore  decay,  and  the  calculated  lower  bound  of  about  0.12  seconds. 
Finally,  the  absence  of  complete  saturation  in  the  viability  even  for  delay  times  as  long  as  50 
seconds,  suggests  that  the  pore  are  relatively  long  lived  in  accordance  with  previous  reports  [44, 45]. 

Numerical  calculations  of  the  cell  deformations  based  on  the  electro-mechanical  stresses 
given  in  the  previous  section  were  performed  to  determine  the  electric  field  effect  on  cellular  shape. 
Results  for  an  initial  spherical  cell  of  thickness  2  nm  having  a  1  pm  radius  (typical  of  E.  coli  cells, 
for  example)  in  response  to  various  electric  field  values  are  given  in  Fig.  7.  Field  magnitudes 
ranging  from  0-70  kV/cm  were  used.  The  Poisson’s  ratio  v  was  taken  to  be  0.2.  The  steady-state 
deformed  cell  shapes  for  positive  z-  and  y-variables  in  the  x=0  plane,  are  shown  in  Fig.  7.  Due  to  •  he 
inherent  symmetry  of  the  problem,  only  the  first  quadrant  is  specified  for  simplicity.  The  shape 
changes  from  a  perfect  circle  for  E  =  0  V/cm,  to  ellipsoidal  with  increasing  eccentricity  at  higher 
fields.  The  corresponding  forces  per  length  N$  (<|>)  and  Ne  (<|>)  are  shown  in  Fig.  8  for  fields  of  20, 50 
and  70  kV/cm.  The  magnitudes  range  from  0  to  about  25  nN/m.  As  reported  in  the  literature  [38],  the 
typical  tension  for  membrane  rupture  is  in  the  range  1-10  nN/m.  Our  results  are  thus  in  very  good 
agreement,  and  show  that  for  applied  electric  fields  of  50  kV/cm  and  higher,  one  can  expect 
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membrane  rupture  simply  based  on  electro-mechanical  considerations.  The  exact  value  will 
obviously  depend  on  the  rigidity  parameter  K  and  the  Poisson’s  ratio  v,  but  the  magnitudes  as 
predicted  by  this  simply  analysis  should  roughly  remain  valid. 

The  deformed  cell  shape,  as  will  now  be  demonstrated,  strongly  depends  on  the  cell 
characteristics.  Changes  in  the  rigidity  parameter  or  the  membrane  thickness  alter  the  force 
distributions,  and  hence,  affect  the  overall  shape.  Calculated  results  of  the  deformed  geometry  for  a 
1  pm  radius  starting  from  an  unstressed  spherical  cell  are  given  in  Fig.  9.  The  membrane  thickness 
ranged  from  2-5  nm  and  various  electric  fields  were  used.  The  curves  of  Fig.  9  show  very  clearly 
that  besides  applied  electric  fields,  the  deformation  is  controlled  by  the  membrane  thickness,  and 
increases  with  “t”.  As  the  thickness  changes  from  2  nm  to  5  nm,  the  geometry  is  modified  from 
spherical  to  ellipsoidal  and  then  begins  to  assume  a  “peanut”  shape.  Based  on  the  trend  evident  in 
Fig.  9,  one  could  qualitatively  predict  an  eventual  shift  towards  a  “dumbbell”  geometry  at  higher 
fields,  or  thicker  membranes,  or  under  conditions  of  smaller  rigidity  parameter,  or  with  a  larger 
Poisson’s  ratio.  Such  calculations  for  strongly  deformed  shapes,  however,  have  not  been  shown  here 
since  the  perturbative  theory  used  in  this  analysis  could  be  called  into  questions  for  large 
deformations. 

Fig.  10  shows  the  bending  moment  M*  (<|>)  and  associated  transverse  force  Q*(<|>)  for  an 
applied  field  of  20  kV/cm  for  an  initial  5  nm  sphere.  As  seen  from  the  curve,  the  magnitude  ofM^, 
(<|>)  is  negligibly  small  and  has  a  nearly  constant  value  of  about  3  x  1  O'12  Newtons.  The  curve  for  Me 
(<)>)  was  nearly  identical  to  that  of  M*  (<t>),  and  so  has  not  been  shown  separately  in  the  figure.  This  Me 
(<j>)  ~  M,),  (<j>)  condition  obtained  here  is  in  keeping  with  a  previous  result  reported  by  Pamplona  and 
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Calladine  [46] .  The  angular  dependence  of  Q^)  from  Fig.  1 0  is  seen  to  be  symmetric  about  <|> = 45°, 
and  also  has  a  relatively  small  value.  Thus,  compared  to  both  N$(<|>)  and  Ne  (<f>),  the  variables  Q$(<J>), 
Me  (<|>),  and  (4»>  can  all  be  neglected  as  has  been  done  in  the  past.  Finally,  Fig.  1 1  shows  the 
fractional  change  in  the  cellular  surface  area  and  volume  as  a  function  of  the  applied  electric  field. 
Two  points  are  evident  from  die  results.  First,  both  curves  exhibit  a  rough  quadratic  behavior.  This 
is  in  keeping  with  some  recent  optical  scattering  experimental  data  [47].  The  exact  magnitudes, 
however,  are  subject  to  the  inaccuracies  and  uncertainty  of  the  material  parameters  such  as  the 
rigidity  K  and  Poisson’s  ratio  v.  Hence,  the  data  of  Fig.  1 1  does  not  lend  itself  to  a  direct 
comparison  with  experimental  data.  However,  the  general  electric-field  dependent  trend  predicted 
here  has  been  shown  to  be  accurate.  A  second  point  about  Fig.  1 1  is  that  the  change  in  cell  volume  is 
larger  than  the  corresponding  areal  variation.  This  is  to  be  expected  as  the  volume  scales  more 
rapidly  than  the  surface  area,  at  least  for  the  simple  ellipsoidal  shapes.  At  higher  electric  fields 
beyond  the  25  kV/cm  value  shown  in  Fig.  1 1,  it  is  conceivable  that  the  areal  variations  become  larger 
as  the  cell  changes  from  an  ellipsoidal  to  a  “peanut-geometry” . 

IV.  SUMMARY  AND  CONCLUSIONS 

A  self-consistent  model  analysis  of  electroporation  in  biological  cells  has  been  carried  out 
based  on  the  coupled  Laplace-Nemst-Planck-Smoluchowski  equations.  The  physical  processes  of 
pore  generation,  drift  and  diffusion  in  r-space  were  all  comprehensively  included.  A  pore-radiiis 
dependent  energy  barrier  to  ionic  transport,  accounted  for  cellular  variations.  This  enables 
predictions  and  qualitative  understanding  of  the  cellular  response  to  short,  high  electric  field  pulses 
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by  taking  account  of  the  growth  and  resealing  dynamics.  The  electroporation  dynamics  in  the 
presence  of  multiple  electric  pulses  and  the  potential  benefits  have  also  been  analyzed. 

It  has  been  shown  that  a  finite  time  delay  exists  in  pore  formation,  and  leads  to  a  transient 
overshoot  of  the  transmembrane  potential  Vmem  •  The  membrane  potential  itself  remains  around  the 
1 .0  Volt  value  that  has  been  reported  in  the  literature.  However,  the  peak  can  exceed  this  value  on  a 
transient  basis.  It  has  also  been  demonstrated  that  pore  resealing  is  a  multi-step  process.  It  consists 
of  an  initial  fast  decay,  followed  by  a  1  O'4  second  delay,  and  then  a  much  slower  pore-closing  with  a 
time  constant  of  about  10'1  second.  This  establishes  a  0.1  millisecond  (ms)  time-window  during 
which  the  pores  are  mostly  open,  and  hence,  the  system  is  most  vulnerable  to  destruction  by  a  second 
electric  pulse.  The  existence  of  such  a  time  window  threshold  for  effective  killing  by  a  second  pulse 
is  amply  supported  by  our  experimental  data  for  E.-coli  cells  involving  two  pulses  with  variable 
delays.  A  sudden  jump  in  the  experimental  cell  viability  has  been  observed  in  going  from  a  0.1  ms 
delay  to  a  1 .0  ms  delay  for  the  second  pulse.  This  is  in  keeping  with  our  simulation  results  that  show 
a  0. 1  ms  time-window  during  which  the  pores  are  mostly  open.  Second,  no  significant  change  in  the 
experimental  viability  was  recorded  in  going  from  a  1 .0  ms  delay  to  a  1  second  delay.  This  again  is 
in  keeping  with  the  prediction  of  slow  pore  decay,  and  the  calculated  lower  bound  of  about  0.12 
seconds.  The  time  constant  for  the  longer  process  also  matches  experimental  results. 

It  has  also  been  shown  that  the  pore  decay  time  for  die  two-pulse  situation  becomes 
significantly  longer  than  a  single  pulse  case.  This  implies  that  applying  two  or  more  pulses  can  slow 
resealing  considerably.  From  a  practical  standpoint,  this  has  applications  for  drug  delivery  over 
prolonged  durations,  or  for  controlled  manipulation  of  the  pore  “open  times”  via  ultrashort,  multiple 
pulses. 
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In  addition,  a  model  analysis  of  cellular  deformation  and  shape  change  in  response  to  an 
applied  electric  field  has  also  been  carried  out.  A  first  principles  approach  based  on  thin-shell  theory 
was  used  to  evaluate  the  deformations.  This  approach  utilizing  the  force  and  moment  equations  has 
the  advantage  that  it  can  potentially  be  extended  to  include  dynamical  analysis  and  temporal 
response.  The  energy-based  schemes  would  fail  in  this  regard.  The  present  calculations 
demonstrated  the  following  features:  (i)  At  low  values  of  the  applied  electric  fields  (as  was 
commonly  the  case  in  the  past),  the  deformed  cells  can  roughly  be  approximated  by  an  ellipsoidal 
shape,  (ii)  However,  for  much  larger  field  magnitudes,  as  have  recently  been  used  [13-16],  the 
deformations  would  be  fairly  significant  and  the  cell  geometry  would  no  longer  be  described 
accurately  by  ellipsoidal  shapes.  For  example,  at  fields  on  the  order  of  50  kV/cm,  a  “peanut-shaped 
geometry  was  shown  to  result,  (iii)  The  results  here  also  demonstrated  that  the  final  shape  depends 
on  the  membrane  thickness.  In  general,  it  was  argued  that  with  decreasing  thickness,  deviations  from 
the  unstressed  shape  would  be  less  severe.  This  has  direct  implications  for  cells,  tissues  and  lipid 
bilayers  in  which  significant  molecular  re-orientation  and  re-structuring  can  occur  upon  the 
application  of  an  electric  field,  (iv)  Values  of  the  surface  forces  obtained  in  the  present  calculations 
were  in  remarkably  good  agreement  with  the  1-10  nN/m  range  for  membrane  rupture  that  has  been 
reported  in  the  literature  [38].  This  lends  validity  and  credence  to  the  present  approach  and 
calculations,  (v)  It  was  also  shown  that,  at  least  for  the  smaller  electric  fields,  both  the  cellular 
surface  area  and  volume  would  change  roughly  in  a  quadratic  manner  with  electric  field,  (vi) 
Finally,  it  was  shown  that  the  bending  moments  are  generally  quite  small  and  can  be  neglected  for  a 
simpler  analysis. 
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Table  I  Parameters  Used  For  The  Theoretical  Model 


PARAMETER 

SOURCE 

VALUE 

D  (m  2  s'1) 

Ref.  19 

5  x  10  14 

Y  (J  m  1 ) 

Ref.  19 

1.8  x  10'11 

r0  ( J  m'2) 

Ref.  19 

10'3 

C  (J 1/4  m) 

Ref.  27 

9.67  x  10 ',5 

Kw(Fm') 

Ref.  19 

80  x  8.85  x  10  ',2 

Km(Fm-‘) 

Ref.  19 

2  x  8.85  x  10  *12 

h(m) 

Ref.  48 

5  x  10'9 

apCFm'2) 

Ref.  27 

6.9  x  10'2 

vc  (m'3  s'1) 

Ref.  36 

2  x  1038 

v«i  (s'1) 

Ref.  48 

10n 

r0(m) 

Ref.  48 

1  x  10'9 

oo  (N  m*1) 

Ref.  48 

5  x  10'2 

oin  (S/m) 

Ref.  15 

0.455 
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o0Ut  (S/m) 

Ref.  15 

5.0 

a  (S/m) 

Ref  15 

1.3 

A  (volt) 

Ref  48 

2.5 

n  =  hi  /h 

Ref  15 

0.15 

iv  (m) 

Ref  34 

65  x  10  '9 

K(J) 

Ref  49 

5  x  10‘2° 
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FIGURE  CAPTIONS 


Fig.  1  Schematic  of  the  model  used  to  represent  a  cell  in  a  suspension  for  potential  calculations. 

Fig.  2  Simulation  result  for  the  pore  density  evolution  with  time  in  response  to  a  4  micro¬ 
second,  10  kV/cm  rectangular  electric  pulse,  (a)  Logarithmic  scale,  and  (b)  Semi- 
logarithmic  graph. 

Calculated  temporal  variations  of  the  transmembrane  potential  in  response  to  a  4  micro¬ 
second,  10  kV/cm  rectangular  electric  pulse. 

Simulation  result  comparing  the  pore  density  evolution  with  one  and  with  two  4  micro¬ 
second,  10  kV/cm  rectangular  electric  pulses.  The  delay  for  the  two-pulse  simulation  was 
taken  to  be  1.2  milliseconds. 

Fig.  5  Simulated  pore  distribution  functions  at  various  times  for  single-  and  dual-pulse  situations. 

Fig.  6  Experimental  results  on  the  viability  of  E.  coli  subjected  to  two  4  ps,  13  kV/cm  electric 
pulses  with  variable  delays.  The  delay  between  the  first  and  second  pulses  ranged  from  5  ps 
to  50  seconds. 

Fig.  7  Calculated  equilibrium  cell  shapes  along  the  y-z  plane  in  response  to  applied  electric  fields 


Fig.  3 


Fig.  4 
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of  0, 20, 50  and  70  kV/cm. 


Fig.  8  Corresponding  forces  per  length  (<|>)  and  Ne  (<j>)  for  applied  fields  of  0,  20,  50  and  70 
kV/cm. 

Fig.  9  Deformed  cell  shape  results  for  various  membrane  thicknesses  and  applied  fields  of  20 
kV/cm  and  50  kV/cm. 

Fig.  10  The  bending  moment  (4>)  and  associated  transverse  force  Q$(<f>)  for  applied  field  of  20 
kV/cm  for  an  initial  5  nm  spheroid. 

Fig.  1 1  Calculated  variations  in  the  cellular  surface  area  and  volume  with  applied  electric  field  for  an 
initial  1  pm  cell  radius  and  5  nm  membrane  thickness. 
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